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The Null Surface Formulation of General Relativity is developed for 2+1 dimensional gravity. The 
geometrical meaning of the metricity condition is analyzed and two approaches to the derivation of 
the field equations are presented. One method makes explicit use of the conformal factor while the 

f—} other only uses conformal information. The resulting set of equations contain the same geometrical 

meaning as the 4-D formulation without the technical complexities of the higher dimensional analog. 

£SJ ■ A canonical family of null surfaces in this formulation, the light cone cuts of null infinity, are 

constructed on asymptotically flat space times and some of their kinematical aspects discussed. A 
particular example, which nevertheless contains most of the generic features is explicitly constructed 
and analyzed, revealing the behavior predicted in the full theory. 
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The Null Surfaces Formulation (NSF) developed by Frittelli, Kozameh, and Newman |l|-|3| gives an alternative 
treatment of Einstein's General Relativity (GR), which emphasizes the dynamical role played in the theory by the 
' null hypersurfaces of the space-time metric. The construction is set up ab initio in terms of generic space-time 
i—t , foliations. These are represented by the level surfaces of a two-parametric function Z(x a ; of the x a coordinates, 
1 on which it is imposed the condition that the foliations become null for some space-time metric g a b(x a ) to be obtained 
as a functional of Z{x a ; £, £). Then, an auxiliary function Q(x a ; £, C) fixes the resulting conformal freedom in such a 
way that the Einstein's field equations be satisfied for g a b(x a ). Such equations are expressed entirely in terms of the 
non-local variables Z(x a ; £, Q and Q(x a ; (, £), with no mention at all of the metric and its associated tensors, which 
O are merely functions of them. 

One important issue in the NSF is how to restrict the freedom of (infinitely many) different families of null hyper- 
$_h . surfaces corresponding to the same metric. For asymptotically flat space-times a canonical family of null foliations 
can be selected by considering past null cones from points at future null infinity, J + . In this case the function Z has 
a second meaning. For fixed values of x a the function u = Z(x a ; parametrically describes the intersection of the 
future light cone from x a with J + . This intersection is called a light cone cut of null infinity and it plays a central 
role in reformulating General Relativity in terms of the free data (representing outgoing gravitational radiation) at 
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It is worth mentioning that the final equations of NSF become a system of partial differential equations for the main 
variables that, in spite of the clear physical and geometric contents of the equations, are technically very involved 
and difficult to analyze. This is best exemplified in an alternative formulation of NSF that decouples the original set 
of equations and produces a conformally invariant equation solely for Z, a desirable feature since Z only contains 
information of the conformal structure [|J . Although the resulting equation is equivalent to the vanishing of the Bach 
tensor, it cannot be explicitly written in terms of Z and its derivatives on a reasonable number of lines. 

Thus, from a purely technical point of view it is desirable to derive an analogous NSF formalism in a lower 
dimensional space. The obvious choice is 3-D gravity since its almost trivial character should simplify the formalism 
and make it easier to investigate all the hoped features that the formulation in four dimensions is required to posses 
in order to be consistent with standard GR. Examples of these are caustics and other singularities of null congruences 
predicted by general theorems, which are at present searched for only in a perturbative approach to the NSF field 
equations 
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Moreover, three-dimensional GR can be used to study many examples of 4-D vacuum space-times with a space 
translation symmetry. By going to the "manifolds of orbits" of the associated Killing field, the problem is reduced 
to 3-D gravity coupled to a scalar field acting as matter source of the 3-D Einstein's equation. Assuming suitable 
fall off conditions on the scalar field, the concept of asymptotic flatness can be stated as close as possible to the 
analogous four dimensional treatment ||. Consequently, for asymptotically flat 3-D space-times, one can introduce 
our canonical family of null surfaces, i.e., the light cone cuts of null infinity. In this way, we extend our formalism to 
4-D space-times which are not asymptotically flat in the usual sense since they include cylindrical or axisymmetric 
waves. 

Finally, from the quantum theory point of view, there is some evidence, at least at a linear level in 4-D Q and 
at a non perturbative for some midi-superspaces in 3-D 0, that in contrast with the quantization in terms of the 
local fields like the metric, the quantization of the non-local variable Z and £1 of NSF, results in operators with 
non-diverging expectation values for semi-classical states. Concrete results have been obtained in 3-D since some 
midi-superspaces can be exactly solved and large quantum gravity effects appear for local fields such as the metric 
||. Recently, it has been shown that the quantum light cone cut admits a semiclassical approximation which is 
stable against small perturbations at transplanckian frequencies [Q. Thus, it appears that this non local variable is 
best suited for semiclassical approximations. To develop a more thorough discussion on these issues one should have 
available the full classical theory. 

This work is organized as follows. In section two we introduce one-parameter families of hypersurfaces in 3-D and 
discuss which conditions should be imposed on these surfaces so that they become null. The main results of this 
section, the metric components and metricity conditions, were originally obtained by E. Cartan and S. Chern many 
years ago |9|,[l0| and recently rederived by M. Tanimoto but we follow an approach which is closely related to 
our work done in four dimensions. In section three we derive the field equations for these surfaces including matter 
sources. As in previous results one can decouple those equations and obtain a conformally invariant equation for the 
null hypersurfaces. In section four we review the notion of asymptotic flatness in 3-D and define the light cone cuts. 
Several kinematical issues which differ from the four dimensional case are pointed out and discussed. A particular 
example of a light cone cut is constructed and analyzed. In the conclusions we outline an approach to handle caustics 
and singularities that come together with null surfaces. 



In this section the NSF formalism, as developed in |2j], is adapted to dimension 2 + 1 and the corresponding 
kinematical structures are constructed. Although the basic framework was outlined in detail in 4-D [Q, it is worth 
repeating the construction in 3-D since it is conceptually equivalent to the higher dimensional case though technically 
less involved and therefore easier to follow. 

We begin with a three-dimensional manifold M , and assume we are given one-parameter functions Z(x a , £) of the 
space-time coordinates x a ; and parameter ( on S 1 running between and 2tt. We also assume that for a fixed value 
of the parameter £, the level surfaces of Z, 



locally foliate the manifold M. 

The statement that the level surfaces of Z be indeed the characteristic or null surfaces of some metric g a b(x a ) for 
M is that for fixed values of x a and arbitrary values of £ they satisfy 



An equivalent geometrical statement is that as the families of foliations intersect at a single but arbitrary point x a we 
demand that the enveloping surface forms the light cone from the point x a . Thus, the parameter £ spans the circle of 
null directions. 

The idea now is to solve eq. (|2j) for the five components of the conformal metric in terms of \7 a Z(x a ^Q. Given 
an arbitrary function Z the problem has no solution since we have an infinite number of algebraic equations (one for 
each value of £) for five unknowns. Therefore we must impose conditions on Z{x a ,C t ) so that a solution exists. The 
solution and conditions are best expressed when written in a canonical coordinate system constructed from knowledge 



II. NSF FORMALISM IN THREE DIMENSIONS. 



Z(x a , C) = const., 



(1) 



g ab (x a )V a Z(x a ,QV b Z(x a ,0=0. 



(2) 




u,u,r) = {Z(x a X),dZ{x a X),d 2 Z{x a ,Q)), 



(3) 
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where <9= denotes the derivative of C holding x a fixed. For each value of C these functions form a coordinate system 
intrinsically adapted to the surfaces. Its associated gradient 8 l a and vector dual basis 8i a are given by 

**-=(<#*)-. e ^(w) a ■ (4) 

We now demand that u = const, is a null surface. Then on that surface oj — const, singles out a null geodesic 
whereas r = const, identifies a point on that geodesic. The next available scalar, d 3 Z(x a , £), determines the conformal 
metric and is constrained to satisfy the metricity conditions. Using eq. (^) and assuming the 6 l coordinate system is 
well behaved, one can obtain x a — x a (9\() and define A(u,cj,r,Q = d 3 Z(x a (0\ (), C)- 

A technical point worth mentioning is that the action of the operator <9, the parameter derivative holding x a fixed, 
acting on a function F(u, u>, r, Q is given by 

d = d c + ud u + rd u + Ad r . (5) 

Thus, the action of the d operator on A is given by (||). 

Likewise, d does not commute with the directional derivatives di=-^. The explicit form of the commutation 
relations, needed later, are 

[d u ,d] = d u Ad r , (6) 
[SU 9} = d u + d w Ad r , (7) 

[d r ,d] = d UJ + d r Ad r . (8) 

The metric components and metricity conditions are obtained by repeatedly operating d on (j^). Taking d on (j^) 
yields 

d{g ab V a ZV b Z) = 2g ab V a ZV b dZ = 2g 01 = 0, (9) 
where we have used that dg ab = 0. If we take again d on the above equation we obtain 

g ab {V a ZV b d 2 Z + V a dZV b dZ) = g 02 + g 11 = 0. (10) 

Repeating this procedure gives 

g ab {V a ZV b d 3 Z + 3V a dZV b d 2 Z) = g 02 d r A + 3g 12 = 0, (11) 

g ab (V a ZV b d 4 Z + AV a dZV b d 3 Z + 3V a d 2 ZV b d 2 Z) = 0, (12) 
g 02 d r dA + 4(<7 n 5 w A + <7 12 <9 r A) + 3g 22 = 0. 

Using the commutator (j^) , together with the already obtained metric components, we can rewrite this last expression 
as, 

g 02 (dd r A - 3d u A - J(a r A) 2 \ + 3g 22 = 0. (13) 

Note that the four nontrivial metric components are proportional to g 02 . It is therefore natural to take g 02 as a 
(■-dependent conformal factor fl 2 and write the metric tensor as g ab = fl 2 h ab = f2 2 /i y 6f8 b with 

/0 1 \ 

= -1 -±d r A . (14) 

\1 -\d r A -\8{d r A) + \{d r A) 2 + d u A J 

The conformal factor Q cannot be an arbitrary function of ( since it is defined as 

n 2 = g ab V a ZV b d 2 Z. 
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Therefore, 

Oil 2 = g ab (V a dZV b d 2 Z + W a ZV b d 3 Z). 

Thus, Q and A must satisfy 

3dfl = fld r A . (15) 

Equation ( |l5| ) is invariant under tt(x, £) - ► f2'(x, C) = f(x)Q(x, C) for an arbitrary /(a;). This freedom is a consequence 
of the conformal invariance of the formulation. 

Note also that nothing prevents us from taking further d derivatives on Since the r.h.s of this equation vanishes 
and all the components of the conformal metric have been explicitly constructed from diA, the next d on (|lj) must 
impose a condition on A so that a conformal metric exists. This equation can be written as 

g ab (5V a dZV b d 4 Z + V a ZV b d b Z + 10V a d 2 ZV ' b d 3 Z) = (16) 
5g n d u dA + 5g 12 d r dA + g 02 d r d 2 A + 10(g 02 d u A + g 12 d^A + g 22 d r A) = 0. 

Using the explicit form of the metric components and the commutation relations and dividing out by the conformal 
factor gives 

M [A] = 2 (d{d r A) - d u A - ^(d r A) 2 J d r A - d 2 (d r A) + 3<9(<9 W A) - 6d u A = 0. (17) 



Conversely, if (|17|) is satisfied, then further d derivatives on (|13j) are identically satisfied. 

We summarize this section as follows: given a function Z(x a , Q such that A((9 J , C) = d 3 Z(x a , C) satisfies M[A] = 0, 
then Z = const, is a null hypersurface of g ab . 



Remark 1: Equation (|T7]) is the main metricity condition, (15) is only used to fix the ^-dependence of the conformal 
factor. In fact, given a function Z(x a ,() such that d 3 Z(x a , "satisfies M[A] = 0, then Z{x a ,(,) = const, is a null 
hypersurface of h ab (x a , (') for an arbitrary value of C, 1 . 

To see this we first show that h ab (x a , (') satisfies 

dh ab = -^d r Ah ab , (18) 
which follows from the definition of h ab . The above equation can be integrated immediately to give 
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h ab (x a X') = exp ^--j^ d r A(x a ,n)dr,j h ab (x a ,Q. 

Thus, if Z(x a X) yields null surfaces for h ab (x a X) it also gives null surfaces for any h ab (x a X') that satisfies (|l§|). 
Since h ab (x a , () is an explicit functional of A we conclude that ( |l7| ) is the only condition to be imposed on Z so that 
a conformal geometry can be constructed on the three dimensional manifold. 

Remark 2: In fact, it is not necessary to assume that g ab = g ab {x) to obtain the explicit form of the conformal 
metric and metricity condition (0). If the metric g ab is such that 

dg ab = »(x,0g ah 

for an arbitrary function u(x, (), we can repeat the steps shown above to obtain the same conformal metric h ab (x a , () 
and metricity condition (|17|). 

Remark 3: Starting from a completely different perspective, E. Cartan || and S. Chern |n| studied the geometry 
of the ordinary differential equation 

d 3 Z = A((,Z,dZ,d 2 Z). 

They found that the condition ( |l7|) has a simple geometrical meaning, namely, that the Cartan normal conformal 
connection defined in the solution space is torsion free jl0|]. It is therefore natural in our case to obtain ([lT]) as a 
metricity condition. 

Remark 4: A very interesting (and open) problem in 4-D is to understand the geometrical meaning of the metricity 
conditions It is natural to conjecture that the complex PDE derived in that case is also the torsion free condition 
for the conformal connection defined in the solution space of the starting second order PDE for Z on the sphere of 
null directions. 
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III. THE FIELD EQUATIONS. 



In this section we derive the field equations for the basic variables of this formalism. In contrast with the 4-D case, 
the vanishing of the trace free part of the Ricci tensor does not contain interesting solutions. Since the Weyl tensor 
vanishes in three dimensions, the solutions of 

Rab — 7^R9ab = 

yield spaces with constant curvature. Thus, we must add a stress energy tensor as a source term on the r.h.s. of the 
above equation. For simplicity we will assume that T ab is a given trace free tensor. If the matter source contains a 
trace part then the trace equation is automatically satisfied by virtue of the vanishing of the divergence of T a b- This 
can be easily seen by writing 

Gab = Q_ ab + -g a bG, 

where G and G ab are the trace and the trace free part of the Einstein tensor respectively. If the trace free equations 

G-ab = £afc> 

are satisfied, then using the vanishing of the divergence of the stress energy tensor and the Einstein tensor one obtains 

V a (G-T)--4V b (G Qb -T ab ) -0. 

We thus start with the Einstein equations and contract the indices a, b with the null vectors 02 a 02 b - Note that 
this operation automatically removes the trace part of the equation. Rewriting the resulting equation in terms of our 
variables yields 

d?n = r rr n. (19) 

At a first glance it appears that this single equation cannot be equivalent to the five components of the Einstein 
equations. Note, however, that (|l^) is true for any value of £. Thus, if we add to ( fl9| ) the metricity conditions which 
gives the change of 17 and A under variations of ( we obtain a consistent set of equations equivalent to the standard 
Einstein equations. The final equations read, 



(E) d 2 r n = T rr n, 

{mi) dVL =-<9 r Afi, 
3 

(mu) Gd u A = -|(9 I .A) 3 + 29 r A9(5 r A) -2d u Ad r A - d 2 (d r A) + Sd^A). t 



(20) 



The above constitute a coupled set of equations for f2 and A . These two functions, however, have completely 
different meaning; whereas A determines the entire conformal structure, only fixes the scale. It is thus desirable 
to decouple the above set of equations and obtain equations that only involve A . By construction, the resulting 
equations would be conformally invariant. 

The approach followed to decouple the field equations is to study the integrability conditions for (mi) and (E) since 
both equations must be satisfied by a single function f2. Therefore, conditions are imposed on A so that a solution 
exists. 

To illustrate the procedure and to discuss the freedom left in the solutions we first consider the case when T rr = 0. 
Using (^]), (0) and (g), we calculate 

[d 2 rl d]Sl = 2d^d r n + 2d r A d 2 r VL + d 2 r A d r VL. (21) 
On the other hand, using (mi) and (E), we have 

[d 2 r ,d]n = d 2 r ( l -d r Ati), (22) 



the r.h.s. of equations (|2l|) and fl22j), gives 



2d u d r tt = \dlAVL- \d 2 r Ad r tl. 
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Applying d r to the above equation, commuting the partial derivatives, and using eq. (E) yields 

d*A = 0. (23) 

Thus, for spaces with constant curvature the field equations for A are, 

d?A = 0, (24) 

6d u A = -^(d r A) 3 + 2d r Ad{d r A) - 2d^Ad r A - d 2 (d r A ) + 3<9(<9 W A ). (25) 
9 

By construction, equation ( pi|) is conformally invariant, ft is thus interesting to ask what tensorial quantity is 
represented by A . It can be easily shown that 

\d A r A = B ab Z a Z\ 
b 

with 

B a b = 2e a " m V[ m S'„] b , 

Sab = Rab — ^Rg a b- 

2V[ m S' Tl ]b known as the Bach tensor in the literature, was first defined and used by E. Cotton |L2j in 1899 to show 
that its vanishing determines the conformal structure of spaces with constant curvature (see Appendix A). 

The solutions to equations (|24|) and ( |25| ) should contain functional degrees of freedom since any null surface com- 
patible with a given Einstein metric should be a solution of those equations. This feature is shown below where we 
discuss the simplest set of equations corresponding to flat space and we obtain the collection of all null hypersurfaces 
in Minkowski space. For this we consider the special class of solutions to the above equations given by A 's that satisfy 



d r A = 0, d^A = -1, 



since this restriction yields a Minkowski metric (see eq. (|14|)). We immediately see that equation (24) is identically 
satisfied whereas the other equation gives 

d u A = 0. 

Thus, the solution can be written as 

A + u = g(Q, 

where g is an arbitrary function on S . Using the relationship between A, u>, and Z we rewrite the above equation as 

d 3 Z + dZ = g((), 

whose solution contains three arbitrary constants. Denoting those constants as (t, x,y) we write the solution as 

Z = < + xcos(C) + ysin(C) + a(C), 

with d 3 a = g(C)- The above solution is the generating function for an arbitrary null hypersurface in flat 3-D space 
p3[ . Its freedom is given by an arbitrary function on the circle of null directions. In the next section we will cut down 
this freedom by finding a canonical family of surfaces. 

Following a similar approach we obtain the field equations for A when matter source is included. Then eq. (^) is 
replaced by (see Appendix) 

dfA = 6d u T rr - 6d r T rt0 + d^A + efT^lW^Q - T ruJ d r VL). (26) 

Note that f2 appears explicitly in the above equation. Thus, it would appear that we are unable to get rid of the 
conformal factor. This is not so. Applying several times the d operator on ( p6| ) gives five independent equations 
involving the three components of V a f2. We can algebraically solve for these three components and leave two extra 
equations for A and the energy-momentum tensor. An explicit illustration of this procedure written in tensorial 
language is given in Appendix B. 
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IV. LIGHT CONE CUTS OF J 



The main results obtained in previous sections, the metricity conditions and the field equations, are valid for any 
family of null hypersurfaces. In this section we introduce a canonical family with a well defined geometrical meaning. 
They represent past null cones from future null infinity. As a necessary step before introducing this family of null 
hypersurfaces we first review the notion of asymptotic structure available in 3-D || . 

A space-time is considered to be asymptotically flat at null infinity if one can attach a boundary with topology 
S 1 x R to the original manifold, a smooth metric field g' ab and scalar field u> on M U J + , such that: 

1. on M u ^ 0, g' ab = cj 2 g ab , 

2. on J+ uj = 0, uj a = V a uj ^ 0, and g' ab uj a uj b = 0, 

3. if uj is chosen such that V a cu a = on J7 + , then the vector field n a = g' ab u>b is complete on J + . 

As in 4-D, the smooth unphysical metric on M U J + is used to discuss the fall off behavior of matter fields near 
J + . By bringing null infinity at a finite distance, one can use the techniques of differential geometry near J + instead 
of taking limits along null directions. Moreover, the space of null geodesies is the same for both the physical and 
unphysical metric, the only difference being in the corresponding affine parametrization of the geodesies. 

There are, however, several differences with 4-D asymptotia. In 3-D the Weyl tensor vanishes and the metric is flat 
outside sources. Thus, the dynamical degrees of freedom are coded in the behavior of matter fields near . 

We now turn our attention to the description of light cone cuts of J + for 3-D asymptotically flat space-times. A 
straightforward method to obtain these cuts is to start with the future light cone from an arbitrary point x a of the 
space time and define a light cone cut of null infinity as the intersection of this cone with the null boundary. Using 
Bondi coordinates (u, £) at J + one can locally describe this intersection as 

u = Z(x a ,()- (27) 

The function Z is called the light cone cut function. It is easy to show that Z has a second meaning. For fixed values 
of u and (, the points x a that satisfy ( ]27j ) lie on the past null cone from the point (u, £) a t J + ■ Thus, the level 
surfaces of Z(x a , form a canonical one-parameter family of null surfaces on the space time. 

In the above construction one makes explicit use of the background metric to first obtain the null cone from x a 
and then intersect this cone with . Since our goal is to regard the cuts as the main variables, we search for a 
construction that does not involve explicit use of the metric. The idea is to define the cuts as the solutions to a third 
order ODE where the source term, A, is obtained from field equations that satisfy the desired asymptotic conditions. 

As was mentioned in the previous section one has available two sets of field equations for A and it is worth reviewing 
both approaches since they involve different variables. The first set of equations, namely eq. (|2(]), are PDE's for 
and A with the matter field T rr as a source term. Imposing regularity and asympotic conditions on the stress-energy 
tensor as discussed in || together with apropriate initial conditions on Q and A one should be able to find solutions 
that preserve the asymptotic structure. For example, f2 and d^A should tend to 1 and -1 respectively as we approach 
to J + . Assuming one has solutions to (|20|) and denoting those solutions as 



fi(it, w, r, C) = Q[T rr ], A(w, u>, r, £) = A[T rr ], 

where the term in bracket denotes functional dependence, one should then use the obtained A to define the light cone 
cuts as the solutions to 

d 3 Z = A(Z,dZ,d 2 Z,() = A[T rr ], (28) 

where the coordinates (u, u, r, £) have been replaced by (Z, dZ, d 2 Z, £). The resulting equation is a third order ODE 
with a 3-dim solution space. Note that the metricity condition is automatically satisfied since it is part of the field 
equations for A . Thus, one can algebraically obtain a metric in the solution space such that the level surfaces 
Z = const, are null surfaces of that metric. 

If we follow the second approach we should start with eqs.(p5|) and ( |26| ) and follow a similar construction to end 
up with (|2§|). The implicit assumption being done here is that one is able to eliminate ft in (|26| ) and obtain a set of 
field equations for A . Although this second approach offers the advantage of having just one variable the first set of 
equations for Q and A appears to be technically more tractable and easier to handle. 

At this point we mention a possible third approach that yields a different equation for the light cone cuts. We first 
derive an analog of the Sachs' theorem in four dimensions (which gives a relationship between the cuts and the shears 



7 



associated with the Bondi and light cone cuts H]). Since the shear vanishes in three dimensions we use the only 
available structure that is left, the divergence of null congruences. To be more precise, we start with the divergence 
p of a null cone with apex at x a . This optical scalar satisfies the following ODE 

%=P 2 + $ oc, (29) 

where s is an affine length and $> 00 = R a bl a l b '. The null cone congruence condition is imposed by assuming that near 

the apex s = s a , p = . Following a similar calculation outlined in Q one can show that at points of intersection 

s s 

between the light cone from x a and the Bondi and light cone divergences are related through 

d 2 Z = PB - pz, (30) 

where ps is the divergence of the Bondi cuts (usually taken to be zero) and pz the solution of (^9|) evaluated at J + . 
Note that pz is a functional of the Ricci tensor <I> 00 . Thus, we can regard ( p9|) and (|3(]) as the field equations for Z. 
Note also that the resulting equation is a second order ODE which in principle is equivalent to the third order ODE 
previously obtained. However, the solution space of this last equation is 2-dim and it is not clear at this point how to 
generate the required extra dimension. A thorough discussion of these issues will be presented in the future. 

A very interesting and difficult problem is to analyze the behaviour of the solutions to the field equations. Since 
these solutions represent characteristic surfaces, they will have caustics and self-intersections. It can be shown that 
both r = d 2 Z and A = d 3 Z diverge at conjugate points. Thus, the field equations as written above are not suited to 
analyze the behaviour of the solutions in the neighborhood of conjugate points. Possible alternative formulations are 
discussed in the last section. 



A. Light cone cuts: an example 



In this subsection we find the light cone cuts associated with particular simple models of space-times. These cuts 
of null infinity exhibit the typical behaviour that one has to deal with in NSF, namely, the presence of caustics and 
self- intersections. 

As we show below, the simplest axisymmetric and static line element, 

ds 2 = -dt 2 + p- 8 M (dp 2 + p 2 dcf> 2 ) (31) 

with a constant M and polar coordinates (t,p,(f>) does not give interesting cuts of J + . Nevertheless it is useful to 
follow the construction since it gives us a hint of how to modify the background space time to obtain non trivial cuts. 

The metric ( |3l|) is interpreted as the space-time of a particle at the origin since its Einstein tensor G a b is given by 
MS(p)V a tV b t fOj. 

After the coordinate change r — ^- with a = 1 — AM, (no relationship with the canonical coordinate r defined in 
previous sections) it can be put in the form 

ds 2 = -dt 2 + dr 2 + a 2 r 2 d(j) 2 , (32) 

which is a locally flat metric with a global deficit angle of 2n (1 — a). Since the solution is not regular at the world- 
line of the particle p = 0, this line must be removed from the otherwise flat space-time. Thus, the integration of the 
geodesic equation is trivial except for the following. For any apex x a there will be one null geodesic that will point to 
the origin and will not reach J + . Therefore, the light-cone cuts at null infinity will be open curves without caustics. 

In order to obtain closed cuts of null infinity, we must then search for a matter source with a non-singular metric. 
The simplest choice is a ring of mass whose metric can be thought of a flat interior matched to a vacuum external 
solution of the form © §: 

-dt 2 + dr 2 +r 2 dcj) 2 , < r < R , 

ds 2 ={ ~ 2 , j-a , / R , .djj ( 33 ) 



-dt 2 + dr 2 + (^^—^j (r + a) 2 d<j) 2 , R < r, 



where the parameter R is the radius of the ring and a is related to the mass. A calculation of its energy momentum 
gives PI, 
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T ^ = ^R(RTaj Sir - R)VatWbt (34) 

Integration on a spacelike surface gives M — —, -, so that M > implies that we must take a > 0. The 

4(i? + a) 

Minkowski space corresponds to a = 0. There is another range of negative a's that also gives a positive mass but 
these paramenters are not continuously connected to a — 0. 

To solve the null geodesies equation and obtain the null cones we introduce the retarded time coordinate u = t — r, 
and rewrite j33| ) as 

ds 2 = -du 2 - 2dudr + x{rfd(j) 2 , (35) 

with 



r, < r < R, 

, (r + a) 
(R + a) 



X (r) = { R (r + a) R<r (36) 



A first integral of the geodesic equations can be obtained by observing that k a k a = and that the scalar product 
between the null vector fc° and the Killing fields \u an d are constant along the null geodesic. 

kaXl = -u-r = -E, (37) 
fc Q q = x(r) 2 4> = L, (38) 
k a k a = -it 2 - 2ur + x{r) 2 4> 2 = 0. (39) 



Then, inserting equations ( |37J) and (|3§|) in (|39j) gives 



2 

x(r) ) 1 



t = ±E\ll- I ) ; (40) 



with 6 = L/E, and from (B7h, (B§) and (M), we get 



2 



i - 



a;(r) 



? = ± , ^ (42) 



dr 



Integrating ( pl| ) and ( ^2|) between r = ro and r — oo, gives the desired light cone cut. The more interesting case 

occurs when ro > R so we will restrict ourselves to this case. We now consider the circle of null directions from the 

point (ito,ro,</>o = 0) and divide this circle in four quadrants. Since the solution has axial symmetry with respect to 

the line that joins ro with the origin we will only consider the upper quadrants which are distinguished by the initial 

value of r(fo)- The integration on the quadrant where r(rp) > is straightforward and does not give any caustics. 

(Along those lines one takes the positive sign in ( |d]) and ( [l2|).) The value of b determines the starting null direction. 

This parameter ranges between 0, when it points away from the origin and b m = x(tq), when r(ro) —0. It is does 

7T . b 

appropriate to introduce a null angle < 6 < — ■ defined by sin 9 = - — . 

Care must be taken in the quadrant where r(ro) < 0. Here b also ranges between b m , and 0, when it points to the 

origin. Then, the range of 9 in this quadrant is given by — < 9 < n. However, there exists a critical value of b, called 

b c (and a corresponding 9 C ) such that for b < b c the geodesies enter the ring. This value of b c — R is obtained from 
the condition that r(R) = 0. 

We start the integration with f < until the turning point r(r t ) — is reached. At this point x(r t ) = b from 
which the value of r t is obtained. Then, for r t < r < oo we choose the positive sign in ( |d]) and ( |4^) to finish the 
integration. 
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The light cone cut of J + from a point with coordinates (uo, r o, <fo = 0) will be given by the following ^-parametrized 
relations between u and <fi: 



u — u = (r + a)(l — cos#), 



(43) 
(44) 



u — u = (r + a)(l — cos#) — 2a\ 1 



r ° + a ■ a 
— sin B 

R + a 



■ = d ( 1 + ^) -2-^arccos(^-±^sin#), 9 C < 9 < ir ; 
V R' i? it + a 



(45) 
(46) 



We now give a brief analysis of the singularity structure of the light cone cuts for this particular case. Without 
involving heavy use of singularity theorems we can see that caustics will occur when 4> is not an injective function of 
6. (If otherwise, we could invert the relation to obtain 6 = 9 {(ft) insert this back into u = U{9{ft)) and trivially show 
that the cut is the graph of a function without self intersections.) 

The graphs <f> vs 6 for fixed values of R = 3 and a = 1 and two values of r = 5, 10 given below serve as illustrations 
of the most general situation. 




FIG. 1. 4> vs 6 for (a) r = 10 and (b) r = 5. 



In Fig. l.a we observe that around 9 = 9 C the function fails to be injective. This case, however, is not analyzed 
in the theory of caustics since at 9 C the function is not differentiable. We thus adopt the criteria that whenever the 
function <j){9) is continuous but not differentiable, we define a caustic point to be such that there is a change of sign 
of the left and right derivatives at this point. It is worth mentioning that the non differentiability around this caustic 
point follows from the fact that the Christoffel symbols, the coefficients of the geodesic equation, have a jump at 
r = R. By smoothing out the shell with a mass density function we can recover differentiability in the solution of the 
geodesic equation. 

For the other value r = 5 (see Fig. l.b), we note the presence of a second caustic point where -377 (#o) = 0. 

d9 

In general, we will always have a caustic point 9 C for any value of R, a, and r and a second caustic point whenever 
r (> R) satisfies the condition 

2a{r + a) < (R + a) 2 . 
We can also plot the corresponding light cone cuts. 
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FIG. 2. Light cone cut for (a) ro = 10 and (b) vq — 5 
du 

As we can see both graphically and analytically, — — is well defined and finite at the caustic points. These caustics 

dip 

are called "cusps". It is rather surprising that, to first order, the singularity structure around 8 C and 8q is similar. 
cP u d? u 

The behaviours of — „ and — , however, are completely different. Whereas at 9q the second and third derivatives 

d<p 2 dep 6 
blow up, at 9 C they are distributional. 



V. CONCLUSIONS 



We close this work with several comments: 



• The cut function Z only gives a local description of a light cone cut. Using the available formalism to describe 
characteristic surfaces [jlTj one can show that a light cone cut will have self intersections (thus Z will not be 
unique) and singular points where d 2 Z and d 3 Z diverge. It is thus desirable to find alternative frameworks to 
describe the cut near singular points and/or globally on J + . One of these alternative approaches is given by 
the method of generating families of Legendre submanifolds Jl7| which might lead to a global description of the 
light cone cuts. 

• Likewise, the null hypersurface Z — const, starts as a smooth surface near J + but later develops conjugate 
points and caustics after the surface enters a region with T a b ^ 0. A similar approach using generating families 
might give a global description of these null surfaces. In particular, there is a need to find a version of the 
metricity condition that will remain smooth around caustic points. 

• We mention again that an interesting and open problem is to study the geometrical meaning of the metricity 
conditions in four dimensions. We conjecture that they are a torsion free condition on the Cartan connection. 

• The most general singularity of a light cone cut in 3-D can be found on the explicit example presented above. In 
addition, we find another singular behaviour that is not consider in the literature since the Legendre submanifold 
is assumed to be at least C 2 . It appears that by relaxing differentiability conditions on the Legendre submanifold 
one should be able to prove generic results for this new type of singularity. 
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APPENDIX A: CONFORMAL EINSTEIN METRICS AND THE BACH EQUATION. 



In this Appendix we derive the field equations for metrics conformally related to Einstein space times. This 
derivation will be restricted to three and four dimensions. Let M be a n-dimcnsional manifold (n=3 or 4) and g ab 
a metric with arbitrary signature. Note that we are not following the same conventions as in the main body of the 
paper where the physical metric is "unhatted" . This is done to simplify the notation for the conformal metrics and 
their related tensors. 

Assume that the metric g ab satisfies the Einstein's equation, i. e. 

Rab ~ -^9ab — Tab 

where T ao is the stress-energy tensor, and R a b, R, the Ricci tensor and the Ricci scalar of g a b respectively. Introducing 

Sab = ^2 [ Rab ~ 2(^1) H ' 



we can rewrite the equations as 



S a b — Sg ab — — T ab . 



1 



We shall express this equation in terms of a conformal metric g a b which is obtained from g ab by rescaling with the 
conformal factor il. (Note that this conformal factor is not the conformal factor for an asymptotically flat spacetime 
as the one given in section IV.) 

g a b = tf'gab or equivalently g ab = fl 2 g ab . (Al) 

We now write the transformation law of various tensor fields under this conformal rescaling. Let us first decompose 
the Riemann curvature tensor of g a b according to 

Rabcd = Cabcd + 2(g a [ c S d ] b ~ gb[ c Sd]a)', 

here we assume that n > 3. In terms of S a b the Bianchi identities can be put as 

V a S ab = V 6 5 (A2) 

and 

V d Cabcd = -2(n - 3)V [a %. (A3) 



Considering the conformal transformation given by (Al), the tensor S a b transforms 



n- 1 

s ab n = s ab n + v a v 6 o - — 5afc v c fw c a (A4) 



Then the conformal Einsteins equation becomes 



rr 1 

s ab n = s ab n + w a v b n - — ga bV c nv c n, (A5) 



where 



If T 

S a b = 77 [ T a b 79ab 

n — 2 \ n — 1 

with T = T ab g ab . Since the stress-energy tensor satisfies V a T a b = 0, for later use, we calculate 

n\/ a T ab = (n- 2)V a QT ab - V h f7T. (A6) 
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We now ask what conditions should be imposed on g a b so that a solution of (A5) exists. If we consider eq. ( A5) as a 
second order differential equation for f2 it is then clea r tha t, for a non trivial solution to exist, integrability conditions 
must be imposed on the metric g a b- Applying V c to (A5), and antisymmetrizing, we get 



v [c s a]b n + v [c ns a]b = v [c s a]b n + v [c ns a]b + v [c v a] v h o - ^Zfl^-iv^v^n). 



Contracting equation (A5) with V a O we find 



and replacing the last term of (A7) by this result, yields 



v [c s a]b n + v [c ns a]b ~ g b[c v d ns a]d = v [c s a]b n + v [c ns a]b + v [c v a] v b n - g b[c v d ns a]d . 

Finally, using the definition of the Riemann tensor, we calculate 

V[ c V a ] Vf,Sl = -C cabd V d Q. + 3fc[ c S' a ] rf V ti il - V[ c £lS a ] b , 

and equation ( |A8| ) becomes 

2(V[ a 56] c n + V[ a f25b] c — g c [ a V d flS b ] d ) — 20V[ a S , b ] c + C abcd Sl d ^l 
Now we consider the case where the stress-energy tensor is pure trace, i.e. 

Tab = —hgab, 

where A is the cosmological constant. This gives in terms of g a b 

AO- 2 



Sab 



(n-2)(n- 1) 



gab- 



Inserting this in the l.h.s. of (A9), we have 



(A7) 



(A8) 



(A9) 



V[ a Si,] c fi + V[ a nS b ] c ~ g c [a^ d ^Sb]a 



AO- 



■jT (-2V[ Q 05fc] c + V[ a Qgb]c ~ 5c[aV 6 ]0) = 0. 



(n-2)(n- 1) 

In dimension n = 3, using the fact that C a bcd = 0, equation ( |A9| ) becomes 

2V [a % = 0. (A10) 
It is easy to prove that the tensor V[ a S b ] c is conformally invariant since 

^[ a S b ] c = V[ a S b ] c tt + V[ f2Sf,] c - g c [ a V d QS b ] d 

= Q(V[ a S b ] c - C c [ a S b ]i) + V[ a flS b ] c - g c [aV d QSb]d 

= ni7 [a S b]c , (All) 

where 

c l ac ^n- 1 (2d{ a v c) n~g ac v l n). 

In order to obtain the integrability condition for n > 3, we use Bianchi identity QA3] ) and equation ( |A9| ) becomes 

1 



n — 3 



V d C a 



6cd 



o, 



(A12) 



where for simplicity in the calculations we have defined w = logO. Applying V a to ( Al2j ) we find 
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V a V d ujC abcd + V a C abcd V d io - -J—V a W d C abcd = 

n — 6 



and using (A5) and (A12), we get 

1 



-R Cahcd H — V Q V Cahcd -V wV a C a b C (i = 0. 



n — 6 n — 6 

Clearly if n = 4 we have the corresponding Bach equation Jl8| 

R ad 

Bbc '■= —y-Cabcd + V°V Cabcd = 0. 



(A13) 



(A14) 



APPENDIX B: THE BACH EQUATION WITH STRESS-ENERGY TENSOR 

In order to handle a more general stress-energy tensor than the pure trace, we may suppose that the stress-energy 
tensor is given and that it d oes not depend on fi. 

We rewrite equation ( A9) in terms of the stress-energy tensor and after using (|A6|) we get 



where 



tabc 



n-2 



V[ a T b ] c - 



2V[ a S'b] c + C a bcdS d ^> — tabc, 

9c[aS d T b ]d V [a Tg h] 



(Bl) 



n-2 



'ujc _ (n - 3) T 

r + (n - 2 )(n-l) V ^' 



We shall concentrate our attention to two cases, when the dimensions are n = 3 and n = 4. In the first case, equation 
(Bl) becomes 



j T-rd ^[a,Tgb]c 

B ca b = 2 I V[ a T b ] c - ,g c [ a V T b ] d h V[ a wT b ], 



Contracting this equation with Tb c and using twice (A6) , we get 



(T bc T bc - T 2 )V a u; = T bc B cab - 2T bc [ X7 [a T b]c - g c[a X7 d T b]d - Zl^lf 



V b T bc T a c + TV C T CI 



Demanding that T bc T bc - T 2 ^ and inserting (|B3j) in (|BJ), we obtain 



(T de T de - T 2 )B cab = 2 (T de T de - T 2 ) ( V [a T b]c - g c[a V d T t 



'-b]d 



T bc T de B ead 



T ac T de B ebd - 2T ae T bc I V[ a T d ] e - g e [ Q V'T d ] 



2T de T ac I V[ b T d ] e - g e [b^7 l T d ]i 



+T bc (V d T de T a e + TV e T ea ) - T ac (V d T de T b e - TV e T eb ) 



When n = 4 equation (Bl) becomes 



CohcdV d W — V d C abcd — t a , 



where 



^abc — V[ a T b ] c 



5c[ a V d T b ] d V[ a Tg b ] c 



2 3 

As before, applying V a and using Einsteins equations we can write 



be? 



V [a wT b]c - — V 



G 



a u 9b]c 



(B2) 



(B3) 



(B4) 
(B5) 

(B6) 
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2-Bab = T a C a bcd — 2 (t ac bV a UJ + y a t a 



be ) 



(B7) 



A simple calculation gives 



2 {t acb V a u; + V°U C ) = 2V Q V [a T b]c - 



9c[aS7 d Tb\d V [Q Tg b]c \ 1 



T 



- ^T ac T\ - -g\ + 



T 



1 3 / 

+ ^T ac R a b - TV c uV b Lo + -V a uV a u iTbc - ^9b. 

+2\7 a Lu\7 a T bc - ^V a LuV a Tg b c + -V 6 wV c T - V a cj\7 c T ab ■ 



Contracting (B5) with C i and using twice (|Aq), we get 

(3C abod C abcd - 6T bc T bc + 2T 2 )V a u - 12V d (C den a T en ) + 3V a (T de T de ) - 6V d (T ae T de ) 

-3V d T de T e a + 5TV d T da - V a T 2 + 2T e a V e T + 

+ 1ZV Cdenl^ a- 



(B8) 



Demanding that 3C abc dC abcd - 6T bc T bc + 2T 2 ^ and inserting ([b|) in (|B7|), we obtain a similar result than for 
dimension 3. 
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